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Polarization properties of turbulent stochastically inhomogeneous ultra- 
relativistic QED plasma are studied. It is shown that the sign of nonlinear 
turbulent Landau damping corresponds to an instability of the spacelike 
modes and, for sufficiently large turbulent fields, to an actual instability of 
a system. 
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1. Introduction 

Working out a quantitative description of the properties of dense strongly 
interacting matter produced in ultrarelativistic heavy ion collisions presents 
one of the most fascinating problems in high energy physics. The simplest 
(albeit not unique) way of putting the experimental data from RHIC [1] 
and LHC [2] into a coherent framework is to describe the essential physics 
of these collisions as a hydrodynamical expansion of primordial quark-gluon 
matter that, after a short transient period, reaches sufficient level of lo- 
cal equilibration allowing the usage of hydrodynamics. The features of the 
experimentally observed energy flow, in particular the presence of a strong 
elliptic flow, suggest early equilibration of the initially produced matter and 
small shear viscosity of the expanding fluid, see e.g. the discussion in [3] 
and [3] devoted to RHIC and LHC results respectively. 

Can stylized features of primordial quark-gluon matter, in particular 
its anomalously low viscosity, be described within a weakly coupled theory, 
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i.e. as a plasma composed of quasiparticles with the quantum numbers of 
quarks and gluons? To address this question let us recall that extensive ex- 
perimental studies of "ordinary" electromagnetic plasma has demonstrated, 
see e.g. [5], that it is practically never observed in the state of textbook 
thermal equilibrium. Realistic description of the properties of experimen- 
tally observed QED plasma is possible only through taking into account the 
presence, in addition to thermal excitations, of randomly excited fields. The 
resulting state was termed turbulent plasma. Collective properties of turbu- 
lent plasmas are markedly different from those of the ordinary equilibrium 
plasmas. In particular, they are characterized by anomalously low shear vis- 
cosity and conductivity, dominant effects of coherent nonlinear structures 
on transport properties. 

Thus it is natural to consider turbulent QCD plasma as a natural candi- 
date for describing the primordial quark-gluon matter in the weak coupling 
regime. Calculation of shear viscosity of turbulent QGP performed in [6] 
has indeed demonstrated that its shear viscosity is anomalously small. In 
the present paper we focus on studying the leading turbulent contributions 
to polarization properties of turbulent relativistic plasma. For simplicity we 
restrict our consideration to the Abelian case. The non-Abelian generaliza- 
tion is briefly described in Section [3l 



2. Turbulent polarization 

2.1. Theoretical formalism 

A weakly turbulent plasma is described as perturbation of an equili- 
brated system of (quasi-)particles by weak turbulent fields F^ u . In the col- 
lisionless Vlasov approximation, the plasma properties are defined by the 
following system of equations (F™, is a regular non-turbulent field): 



d^-eq {F*+F%,) 



_d_ 



fip,x,q) = 



0" (F* + F%) = j u (x) =e^2 [dppuq f(p, x, q). 

q,s J 



(1) 



The stochastic ensemble of turbulent fields is assumed to be Gaussian and 
characterized by the following correlators: 

<i£) = 0, (F T ^(x)F T ^'(y)) = K^'-»'(x,y) (2) 

In the present study we use the following parametrization of the two-point 
correlator K^ u ^ v (x,y) [6]: 
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Turbulent polarization arises as a (linear) response to a regular pertur- 
bation that depends on turbulent fields. It is fully characterized by the 
polarization tensor H lJ ' u {k) defined as a variational derivative of the aver- 
aged induced current (j M (/c \F r ,F t )) f t over the regular gauge potential 
A R : 



n-W = gfflg (4) 

(f(k\F R ,F T )) FT =eY, ( dP Pu q(Sf(P,k,q\F R ,F T )) FT (5) 

q,s J 

Let us rewrite the kinetic equation in (pQ) in the following condensed form 



f = r + Gp^F^f, G= n „^\ - T , (6) 



eq 

i((pk) + le) 



where f eq is a distribution function characterizing the original non-turbulent 
plasma and introduce the following systematic expansion in the turbulent 
and regular fields: 

*/=EEA^/mn, i^=££p m r n «, (7) 

m=0 n=0 m=0 rt=0 

where powers of p count those of F R and powers of r count those of -F r . 
Turbulent polarization is described by contributions of the first order in 
the regular and the second in the turbulent fields. The lowest nontrivial 
contribution to the induced current © is thus given by 5fi%. We have 



where 



5/ ^ 5/htl + (<5/i2)i + (5/12)11 
5/htl = Gp^d^ p r 

(5/i2)i = Gp^F^d v>p G P ^F&'d„ />p Gp p Fg}d a!P r* 
(5/i2)ii = Gp^F^d u>p G Pfl ,F^'d u/jP Gp p F^)d a , p n 

2.2. Nonlinear Landau damping and instabilities 

Generic expression for the polarization tensor taking into account tur- 
bulent effects can be written as 

n,>,k;0= r<5«-^W(a;,|k| ; I) + ^H L (u, |k| ;l) (8) 
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where I = \f2{ra) / y 't 2 + a 2 . Both longitudinal and transverse components 
can be presented as a sum of Hard Thermal Loops (HTL) contributions and 
the gradient expansion in the turbulent scale I: 

!irL ( W ,k)+n^(u;,k| I) 



n L(T) (w,k ; z) = n 
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x = lj/ |k| and the standard HTL contribution 
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(10) 



The computation of turbulent polarization was carried out to second order in 
the gradient expansion [7\ . In what follows we restrict ourselves to discussing 
the leading contribution to the imaginary part of the polarization function 
corresponding to the turbulent modification of Landau damping in (|10|) : 
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The functions Im <^l T aim ^ m Xlt are shown in Fig. Q] The conclusions 
following from Fig. [T] can be formulated as follows: 

1. Timelike domain. From Fig. [I] we see that the sign of the imagi- 
nary part of the turbulent contribution to the polarization operator 
in the timelike domain x > 1 is negative and corresponds to turbu- 
lent damping of timelike collective excitations. This refers to both 
transverse and longitudinal modes. As the HTL contribution in this 
domain is absent, this turbulent damping is a universal phenomenon 
present for all u, k such that uj > k and all values of the parameters 
involved (I, (B 2 ), (E 2 )). The turbulent damping leads to an attenua- 
tion of the propagation of collective excitations at some characteristic 
distance. 
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Fig. 1. The polarization functions (67r 3 / 2 /e 4 ) Im [</> ^ (x)] (solid lines) and 
(67r 3 / 2 /e 4 ) Im [x ^^(a;)] (dashed lines). Left: transverse response; right: longi- 
tudinal response. 

2. Spacelike domain. The situation in the spacelike domain x < 1 is 
more diverse. In contrast with the timelike domain the gradient ex- 
pansion for the imaginary part of the polarization tensor starts from 
the negative HTL contribution corresponding to Landau damping. As 
seen from Fig. [T]the imaginary parts of turbulent contributions to the 
longitudinal polarization tensor are negative and are thus amplifying 
the Landau damping. The most interesting contributions come from 
the turbulent contributions to the transverse polarization tensor. We 
see that the electric contribution Imf^ 1 ^)] in the spacelike domain 

is positive at all x while the magnetic contribution Imfx^^)] is neg- 
ative for x < x* ~ 0.43 and positive for x > x* . This means that the 
turbulent plasma becomes unstable for sufficiently strong turbulent 
fields. 

3. Comments and conclusions 

Let us briefly discuss some relevant points: 

1. The above presented results are obtained in the framework of a pertur- 
bative expansion based on two crucial assumptions. First, one assumes slow 
temporal evolution of the distribution function due to particle interaction 
with turbulent fields thus neglecting the corresponding fo n contributions. 
Second, changes in the distribution function are treated as small. This, in 
turn, means that turbulent fields should be small enough. In this sense the 
reliable results refer to small modifications of Landau damping but, as the 
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onset of turbulent instability takes place for parametrically large fields, this 
result should be considered as a qualitative indication. 

2. The observed instability can be termed "secondary" because the turbu- 
lent fields themselves result from some "first level" instabilities. The origin 
of the effect is in turbulent stochastic inhomogeneity and thus similar to 
stochastic transition radiation, which vanishes in the limit I — )■ 0) [8]. 

3. The non-Abelian generalization of the above-described results for the 
imaginary part of the polarization tensor leads to identical expressions, the 
only difference being in trivial color factors, just as in the HTL case. 
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